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SUMMARY 


An analytical model of an Annular Momentum Control Device (AMCD) 
laboratory test model magnetic bearing actuator with permanent magnet flux- 
biasing is presented. An AMCD consists of a spinning annular rim which is 
suspended by noncontacting magnetic bearings and driven by a noncontacting 
linear electromagnetic spin motor. The actuator is treated as a lumped- 
parameter electromechanical system in the development of the model. 


INTRODUCTION 


The basic concept of the Annular Momentum Control Device (AMCD) is that 
of a rotating annular rim suspended by noncontacting magnetic bearings and 
powered by a noncontacting linear electromagnetic motor. A detailed discussion 
of the rationale for the AMCD configuration and its potential applications are 
presented in reference 1. Earth-based energy storage applications of the 
concept are discussed in references 2 through 4. 

In order to investigate any potential problems in implementing the AMCD 
concept, a laboratory test model AMCD was designed and fabricated under 
contract. The laboratory model has been delivered and preliminary tests have 
been performed. Reference 5 presents a detailed description of the laboratory 
model and reference 6 presents results of static and low-speed dynamic tests 
which include spin motor torque characteristics and spin motor and magnetic 
bearing drag losses. Reference 6 also briefly discusses permanent magnet 
flux-biasing. Previous analyses of the laboratory model AMCD (including refs. 

5 and 6) have utilized a linear model of the magnetic bearing actuator. One 
limitation of the linear model is that it is accurate only in a small region 
about a given operating point in the bearing gap. The purpose of this report 
is to develop a nonlinear analytical model of the AMCD laboratory model 
magnetic bearing actuators which can be used .in computer simulations to 
investigate the effects of large rim' motions in the bearing gaps. 
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cross sectional area of magnetic bearing element core 
magnetic flux density , 

magnetic flux density in permanent magnet element 

residual magnetic flux density in permanent magnet 
element 
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electric field intensity 

electric field intensity measured with respect to a 
contour that may be moving 

net force produced by magnetic bearing actuator 

force produced by lower magnetic bearing element 

total force acting on suspended mass 

force produced by upper magnet bearing element 

force of electric origin 
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lower magnetic bearing element gap 

upper magnetic bearing element gap 

magnetic field intensity 

magnetic field intensity in permanent magnet element 

defined by equation (29) 
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position gain 
electromagnet gain 

equivalent permanent magnet stiffness 
rate gain 

defined by equation (39) 

defined by equation (39) 

length of permanent magnet element 

magnetization density 

mass 

number of turns in electromagnet winding 

normal vector 

laplace variable 

voltage 

weight 

total energy stored in magnetic field 
coenergy 

magnetic susceptability 
displacement 

bearing gap with rim centered 

position sensor output 
defined by equation (41) 
flux linkage 
permeability 



permeability of permanent magnet material 
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A 


permeability of free space 


magnetic flux 
equal by definition 


SINGLE ELEMENT MODEL 

In this section the force equations of a single magnetic bearing element 
are developed. The element is treated as a lumped-parameter electromechanical 
system and the magnetic field equations used in the development are presented 
in Appendix A. For a detailed treatment of lumped-parameter electromechanical 
systems see reference 7. 


Electromechanical Coupling 

Electromechanical coupling is modeled as a lossless magnetic field system 
which couples an electrical terminal pair with variables i and X and a 
mechanical terminal pair composed of a node x acted on by an electrical 
force fg. Flux linkage, X, is discussed in Appendix B. All mechanical 

energy storage and electrically dissipative elements are assumed to be 
connected to the terminals externally leaving only elements that store energy 
in the magnetic field. Because of this the assumption is made 
that the force of electric origin, f , has a single value for a given i 
and X. That is. 


f = f (i,x) 
e e 


( 1 ) 


Energy Considerations 


Let be the total energy stored in the magnetic field. The 

conservation of power can be written as 


iv = 


'e dt 


( 2 ) 


where iv. is the time rate of change of energy (power) going into the 

dx 

electrical terminals and ^ is the resulting time rate of change of 

energy going into the mechanical terminals. Or, in terms of energy, equation 
(2) can be written as 


dW 


m 


4 


dt 


. dX 
^ dt 


f ^ 

e dt 


(3) 



where the relationship of equation (B5) has been used for v. Eliminating 
dt results in the conservation of energy 

dW = i dA - f dx (4) 

Next, choose the independent variables A and x which determine i and 
f^ through the terminal relations. The change in due to changes in A 

and X can then be found by integration of the above equation. This will be 
a line integration in two dimensional variable space and since the system is 
conservative, the stored energy is a function of its state and is independent 
of the path taken to reach that state. Choosing path 1, 2, 3 in figure 1 
results in 



S x^ f \ 

+ 1 i(A,x^)dA 


(5) 


By using this path and starting with the electrical terminal values zero 

instead of an arbitrary value, W can be determined from the electrical 

m 

terminal relations only. That is, ivhen A^ = 0 there is no force of 
electrical origin (f^(0,x) = 0 in the integration). 

Force-energy relations .- From the section above it is clear that can be 

expressed as a function of the independent variables A and x 

= WJX.x) (6) 

Taking the total derivative- 

9W 9W 

dW^ = ^ dA + J dx (7) 

m oA ox 
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and subtracting this equation from the previous equations for dW results 
in 





0 


( 8 ) 


Since X 
values . 
Thus, 


and X are independent variables, dX and dx can have arbitrary 
This means the coefficients must be zero to satisfy the equation. 
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W CX,x) 
m 

9X 


(9) 


and 


f 

e 


9W (X,x) 
m 


3x 


( 10 ) 


Force coenergy relations .- In the present development it is more convenient 
to use i instead of X for an independent variable. From above 


dW = i dX - f dx 
m e 

X = X(i,x) 


and 


f = f (i,x) 
e e 

By defining the coenergy to be 

W Axi - W 
m m 


then 


which simplifies to 


dW = d(Xi) - idX + f dx 
m e 


dW* = Xdi + f dx 
m e 


Since W» = W‘ (i,x) 
m m ^ 


3W» 3W 

dW = di + dx 

m 3i 9x . 


Subtracting equations (13) from (12) results in 

/ 3W \ / 9W \ 


= 0 


( 11 ) 


( 12 ) 


(13) 


(14) 
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from which is obtained the force of electrical origin as a function of current 
and displacement 

3W| (i,x) 


m 


3x 


(15) 


Simplified Magnetic Bearing Element 
with no Permanent Magnet 

Consider the magnetic bearing element shown in figure 2. 



Assumptions : 

(1) 

Magnetic material has infinite 


(2) 

Negligible fringing 


(3) 

Negligible leakage 

From equation 

(A7) 



B = 

yH 


For a finite B inside the magnetic material, with y H 0. This means 
that all the magnetic field intensity is across the air gaps. Using equation 
A1 and the contour shown 



d£ = H^g + H2g = 



By symmetry it can be seen that 
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2g 


(16) 


(17) 
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The flux density through the core is the flux through one of the gaps and is, 
from equation A7, 


y Ni 

B = y H = y H. = — 

o 1 2 2g 


(18) 


The flux (J) in the core becomes 


(f) = BA = 


y^Ni A 
2g 


(19) 


where A is the cross sectional area of the core. This flux links the N 
turn winding N times so that from the definition of X (equation B4) 


X = N 




da = N (}) 


( 20 ) 


where s is the surface enclosed by the winding. Then 

y AN^i 
o 


2g 


( 21 ) 


From the previous definition of coenergy (equation (12)) 


dW» = Xdi + f dx 
m e 


or between the points (i^,x ) and (i, ,x, ) 

3-3. DO 




Integrating from i = 0 and substituting g for x results in 

.i 




w; = \ x(i,g)di 


since 


fg(o,g)dg = 0 


The bar represents the inte^gration variable and (i,g) represents the end 
point of the integration. 


W then is 
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m \ 2g 
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(23) 


8 



From the force-coenergy relations (equation (15)) 


or 
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3W;(i,g) 
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(y^A(Ni)' 




y^A(Ni)' 


4g 


( 24 ) 


Simplified Magnetic Bearing Element 
with Permanent Magnet 

Consider the magnetic bearing element with a permanent magnet in the core 
as shown in figure 3. 



Magnetic Material 
Permanent Magnet 


Figure 3.- Magnetic bearing element with permanent magnet. 
A good approximate model of a permanent magnet is (from reference 7) 


B = y H + B 
m mm o 

Taking the line integral around the contour 




da 


From equation (25) 


m 



(25) 

shown results in 

+ H-g + H ^ = Ni (26) 
2 m 

(27) 
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Equation (26) becomes 


B 


B 


Hig ^ H^g = Ni 

m m 


The term — i is constant and can be represented by 
^m 

magnetomotive force (mmf) 

B A 

NI 

% ° 

As before, = H and B in the air gap = B 

Then, 


B = y H = B 
^o m 


and 


H becomes 


(y H)£ 

H(2g) + = N(i + I^) 


H = 


N(i + I ) 
o 

2g + — £ 

% 


(28) 

an equivalent 

(29) 

in the core. 

(30) 

(31) 

(32) 


and B becomes 


B = y. 


N(i + I^) 
y 

2g + — £ 


(33) 


X can then be determined as 


Next, define 


X = N(J) = NBA = 


i" = i + I 


y AN^(i + I ) 
0 o 

y 

2g + — S, 

^m 


Integrating from- i'’ = 0 and substituting g for x 
of ' 




m 


(34) 


(35) 

results in a 

y AN^i"^ 
o 
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coenergy 

(36) 
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then becomes from equation (15) 
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y AN^(i + I 
o ^ o 



( 37 ) 


SINGLE ACTUATOR MODEL 

A single AMCD axial bearing actuator, for control along a single axis, 
consists of a pair of magnetic bearing elements, with permanent magnets, 
mounted in the configuration shown in figure 4. 


The bearing elements are connected in a differential configuration. 

That is, for a given input the amplifier driver shown in the figure produces 
current in a direction to aid the permanent-magnet-produced flux in one 
element while at the same time- producing equal current in a direction to 
subtract from the permanent-magnet-produced flux in the other element. This 
results in a net force produced on the suspended mass in a direction dependent 
on the polarity of the input to the amplifier driver. If the force produced 
by the upper element is defined as and that produced by the lower element 

as and up is taken as the positive direction, the net force produced by 


an actuator as a function of current and gap becomes 
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F = p AN 
L ^o 


y. 


[I + i]‘ 
o 


[I - i]‘ 
o 


\ 


y 12 
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(38) 


The multiplication of — S- by 2 results from the fact that there are two 

permanent magnet wafers in each AMCD magnetic bearing element (reference 5) . 

By making the following definitions, equation (36) can be written in a simpler 
form: 

y AN^ = K, 
o 1 



(39) 


where x^ is the bearing gap with the rim centered and x is the displacement 
of the rim with respect to x^ (up is positive displacement). 
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Equation (38) then becomes 






[I + if 

, - X + Kj]' 


[I - i]^ 


[\ ■+ X + K^]' 


Finally, define 


- A 

= X + K- 
0 o 2 

The equation for force then becomes 
K 



[I - if 

r- i2 

[Xo - X] 


(40) 


(41) 


(42) 


This equation can be linearized for small deviations about a given operating 
point by using the combination of linear gains 


F = Kgi 

U IX A 3F 
where and 

B 

(reference 8) . Using 


K X 
m 


(43) 


Km = ^ and are evaluated at the operating point 

equation (42) these gains become 


and 


= "l 


I x^ + 2x xi + I X ^ 
o .0 o o 

- 2 2 
X - X 
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K / [I + i]^ [I - i]^' 
K = ^ ^ . 

m 2 \ r“ r*“ 


[Xo - X]- 


[Xo - X]- 


(44) 


(45) 


In the case of the AMCD laboratory model, the parameters and can be 

obtained from hardware data. At an operating point of .635 mm (.025 inches) 
above center, K^ is 13.79 N/amp (3.1 Ib/amp) , and at i = 0 is 

59.54 N/mm (340 Ib/inch) (reference 5). x^ is the gap with the rim centered, 

2.54 mm (.10 inch), plus the thickness of one of the permanent magnet wafers, 

1.02 mm (.04 inches) (reference 5), since — ~ 1. This leads to three 

equations "((42), (44), and (45)) and three unknowns (I^, K^, and i at the 

nominal operating point) . F in equation (42) is one-third of the rim weight 

(16.5 lb). Solution of the equations yields I = 14.97 amp, K = 9.9 N-mm^/ 

2 -3 2 2 ° ^ 

amp (3.45 (10) lb-in /amp ), and i at the operating point = 2.99 amp. 

Solving equation (45) using i at the operating point yields K = 73.67 N/mm 

(420.7 Ib/in) . Equation (42) then becomes 
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F = 2.475 


C46) 


/ [14.97 + i1^ _ [14.97 - i1^ \ 

\[ 3.56 - x]^ [ 3.56 + xf) 

and (42) linearized around an operating point .635 mm (.025 in) above the gap 
center becomes 

F = 13.79 i + 73.67 x ( 47 ) 

With a finite position gain, the equilibrium position of the rim will be 
below the sensor zero. A block diagram representation of a single bearing 
control loop is shown in figure 5. 


1 



Figure 5.- Single bearing control loop. 


At equilibrium with the position command, x , set to zero, the total 
force, F^, from the figure becomes ^ 

p . 2.475 - i)^ \. „ , 0 

\ [ 3.56 - x]^ [ 3.56 + x] y 

where W represents one-third of the rim weight. The current becomes 

i = -Ka ^ = Ka'(Xo - X) (49) 

Using the values of x^ and W from above and a value of 19.57 amp/mm 
(497.1 amp/in.) for (reference 5), equations (48) and (49) can be solved 

simultaneously to yield an equilibrium position of x = .427 mm (.0168 in.) 
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and an equilibrium current of i = 4.08 amp. Substituting these values into 
equations (44) and (45) yields Kg = 13.03 N/amp (2.93 Ib/amp) and 

= 67.88 N/m (387.6 Ib/in) . Equation (42) linearized around the equilibrium 

point becomes 

F = 13.03 i + 67.88 x (50) 


Concluding Remarks 

This paper has developed a nonlinear analytical model of an annular 
momentum control device (AMCD) laboratory test model magnetic bearing 
actuator with permanent magnet flux-biasing. The actuator is treated as a 
lumped-parameter electromechanical system in the development. The model can 
be used in computer simulations to investigate the effects of large AMCD rim 
motions in the magnetic bearing actuator gaps. 
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APPENDIX A 


Magnetic Field Equations 

This appendix presents the magnetic field equations used in the 
development of the analytical model of the AMCD magnetic bearing actuator 
presented in the main text of this report. For a development of these 
equations, see reference 7, Appendix B. 


The magnetic field equations used are: 


H . dK = ^ 



• nda 



nda 


(Al) 

(A2) 

(A3) 

(A4) 


where 


0c ^ 


Pq (H + M) 


(A5) 


is the line integral around a closed contour 


integral over a surface enclosed by a contour 


•0s ' 


■ 0 ^ 


is the surface 


is the surface integral 


over a surface enclosing a volume, n is a normal vector, H 


field intensity, B is magnetic flux density, is current density. 


is magnetic 

E' 


is electric field intensity measured with respect to a contour that may be 
moving, is the permeability of free space, and M is the magnetization 

density which accounts for the effects of magnetizable materials. M is 
usually defined as 


M = X H 
m 

where X^ is magnetic susceptibility. If the permeability of a material is 
defined as 

y = (1 + X^) (A6) 

then equation (A5) can be written as 


B = pH 


(A7) 
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Appendix B 


Flux Linkage 


In this appendix the concept of flux linkage is discussed. For a more 
detailed discussion, see reference 7. 

Consider a perfect conductor with terminals a and b as shown in 
figure B1 where the contour between a and b outside the conductor encloses 
no changing magnetic flux. 






No Changing 
Magnetic Flux 


Figure Bl.- Perfect conductor with terminal voltage v. 


The terminal voltage v can be written as 

a 

V = - ^ E • d£ (El) 

b ' ' 

where the line integral is taken from b to a outside the conductor, 
by definition a perfect conductor can support no electric field 


Since 




dZ = 


(B2) 


18 



From equation (A4) in Appendix A 


V 




nda 


Define flux linkage X of the conductor as 




nda 


V then becomes 


V 


dA 

dt 


(B3) 


(B4) 


(B5) 


From equations (Al] through (A4) in Appendix A, it can be seen that flux 
density B will be a function of terminal current (equation (A3) is used to 
relate terminal current to current density in the system and then (Al), (A2), 
and (A5) used to solve for B), material properties (equation (AS)), and 
system geometry. Therefore, A is also a function of terminal current, 
material properties, and system geometry. If the assumption is made that 
geometry is fixed, except for one movable part whose position can be described 
instantaneously by a displacement x with respect to a fixed reference, and 
that M is a function of field quantities alone then 


and 


A = A(i,x) 

dt 9i dt 9x dt 


(B6) 

(B7) 


where terminal voltage is shown to be a function of rate of change of current 
in the circuit and relative motion of the circuit. 
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